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ABSTRACT 

We have obtained the first large sample of accurate temperatures for 
clusters at z > 0.14 from ASCA. We compare the luminosity temperature (L-T) 
distribution for these clusters with the low redshift sample of David et al (1993) 
and find that there is no evidence for evolution. We also find that the intrinsic 
variance in this relation is roughly constant with redshift. Additionally there is 
no detectable change in the relationship of optical velocity dispersion to X-ray 
temperature with redshift. Most cosmo logical simulations driven primarily by 
gravity predict substantial changes in the L-T relation due to the recent rapid 
growth of clusters. Our results are consistent either with models in which the 
cluster core entropy is dominated by pre-heating, or with low Q models in which 
cluster structure does not evolve strongly with time. The intrinsic variance in 
the L-T relation at a fixed redshift can be due a variety of possibilites e.g. a 
change in the baryonic fraction from cluster to cluster, variation in the fraction 
of the total energy in the system arising from shocks or supernova heating or 
variations in the emission measure distributions in multiphase gas. 

Subject headings: galaxies:clusters:general - X-rays:galaxies - 
cosmology:observations 
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1. Introduction 

Clusters of galaxies are the largest relaxed systems in the universe and as such provide 
a strong test of theories of the origin of and evolution of structure. A critical measurement 
( [Frenk et al 19901) is the relationship of the mass of a system to its temperature. This 
can be best probed by comparing the evolution of the X-ray luminosity (which is related 
to the baryonic mass) to the x-ray temperature, which is roughly proportional to M 2 / 3 . 
Calculations of the evolution of these quantities (of., |Kang et al 1994] , |Kaiser 199 1| ) indicate 



that, in a closed universe, where gravity dominates the evolution of the gas and the dark 
matter, the relationship between luminosity and temperature should evolve strongly at all 
redshifts, with significant changes occurring at z < 0.5. However the nature and amplitude 
of the evolution is dependent on the details of the cosmological model ( Pryan et al 1994]) . 

For example, evolution is less in the CDM+ A=l model than in a pure CDM universe 
and is also less in a low density universe, or one in which the entropy of the gas is relatively 
large at high redshift ([Kaiser 1991| , [Bower 1997| , [Navarro, Frenk and White 1995| ). This 



evolution is not necessarily monotonic. Cen & Ostriker (1994) show that in a CDM+A=1 
model there is an almost constant mean cluster temperature from z ~ 0.3 to the present 
while in a CDM model there is a sharp decrease in the temperature from z = 0.0 to z = 0.3. 



At low redshift the L-T relation is well measured ( Mushotzky 1984] , [Edge fc Stewart 



[199TI |David et alT993|) . While the correlation is good, there is a fair degree of scatter 



( |Fabian et al 1994]) , with cooling flow clusters being rather more luminous for a given 
temperature. The theoretical prediction for the relationship between luminosity and 
temperature in a CDM universe dominated by gravity ( [Frenk et al 1990] , |Kang et al 199 4 



Fvrard 1990| ) is somewhat flatter than observed. There are various possible explanations 



for this such as systematic variations in the baryonic fraction with cluster mass (Pavid et] 
al 1993j) or pre-heating ([Fvrard 1990|). 
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At low redshifts the relationship between the cluster velocity dispersion (ay) and the 
gas temperature (|Bird, Mushotzky fc Metzler 1995| , Pachall fc Lubin 1994| , |Girardi et al 



1996|) is very tight, consistent with both tracing the same potential. This strong agreement 



is not generally predicted to continue at higher redshift. 

ASCA data allow, for the first time, a measurement of the L-T and the oy-T 
relationships at z £j 0.1 from a large sample of clusters. Thereby placing strong constraints 
on theories of cluster evolution. In this paper we use a sample of 38 clusters of galaxies from 
z ~ 0.14 — 0.55 obtained from the ASCA archives. This sample, while not homogeneously 
chosen, is large enough to constrain this relationship at (z) ~ 0.3. 



The only previous attempt to measure the L-T relationship at z £ 0.1 ([Henry, Jiao 



& Gioia 19941) was consistent with no evolution, but with rather large errors. However, 



this study relied on summing a large number of low signal to noise data sets from the flux 
limited Einstein medium survey and was therefore fundamentally different in character. 
The ASCA data have well determined temperatures for each cluster but are drawn from a 
heterogeneous sample. 



2. Observations 

We have extracted the spectra from 38 clusters of galaxies at z > 0.14 from the ASCA 
archive (Table 1). The sample simply consists of all clusters available before Nov. 1996 
for which the data are publicly available, the temperature uncertainty is 5T/T £ 0.3, and 
there is no strong substructure visible in the ASCA image. In our analysis we allowed the 
cluster abundance and galactic absorption to be free parameters. For clusters at z ~ 0.3 
this effectively requires cluster flux x ~ 4 x 10 -8 ergs cm" 2 . 

Examination of the data (Figure 1) shows a relatively strong selection bias against 
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clusters of L\> \ < 10 45 erg s _1 at z ~ 0.3. This comes from the desire of many of the 
proposers to obtain higher quality spectra, requiring bright sources and/or longer exposures 
than normal with ASCA, combined with the relatively small number of known clusters in 
this luminosity/redshift range. 

At z ~ 0.3 clusters of < 10 45 are below the threshold of the ROSAT all sky 
survey and mainly come from the EMSS data base QGioia et al 1990|) and serendipitous 



ROSAT sources. At flux levels £ 10 12 erg s x cm 2 the ASCA data analysis is relatively 
straightforward (see [Mushotzky fc Loewenstein 19971 f° r a discussion of the higher signal 



to noise objects in this sample). We have used the latest calibrations including the gain 
change in the SIS detectors. 

Quoted errors are 90% confidence for one parameter (x 2 +2.71). We use go — and 
Hq = 50 km s _1 Mpc -1 in this paper. Our derived temperatures agree very well with 
those obtained by |Tsuru 19961 , f° r common objects and the previously published results 



of ISchindler et al 1997[ |Allen et al 19961 , |Bautz et al T99l and ponahue 199^ . 



We have integrated the spectra over 3-6 arcmin in radius, depending on cluster redshift, 
in an attempt to get the average spectrum in a fashion similar to those obtained by 
non-imaging proportional counters for lower redshift objects. Thus our values are directly 
comparable to those of the non-imaging experiments. The effects of cooling flows ( |Allen et| 
|al 1996|) , mergers ( [Henriksen fc Markevitch 1996|) or non-isothermal profiles ( [Markevitch 



1996|) can change the mean temperature by up to ±20%. For example, in A1835 which 



has a very large cooling flow, the mean kT=8.15 keV (Table 1 and |Allcm Fabian fc Kneib 



1996|) changes to 9.5 keV with the inclusion of a cooling flow in the spectral fit. The 
uncertainties in the temperature for these high z systems are similar to or less than the 
errors obtained from EXOSAT ( |Edge fc Stewart 1991| ), HEAO-1 flMushotzky T984[ ) or the 



Einstein MPC ( pavid et al 1993Q proportional counters, and thus form a good comparison 
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sample. We have added a few low z clusters to the David et al (1993) compilation (A3158 
and A3581) and obtained accurate ASCA temperatures from our own analysis and the 
literature for those low z clusters with large temperature uncertainties in the David et al 

(1993) compilation (A399, A401, A3112 and A3391) to obtain an essentially flux limited 
low z sample resulting in a total sample of 102 clusters from z ~ 0.01 to z ~ 0.55. 

3. Analysis 

Comparison of this sample (z > 0.14) with the enhanced David et al (1993) sample 
of low redshift clusters (z < 0.1) shows no evidence for a change in the L-T relationship 
(Figure 1) over the whole luminosity range covered (log 10 L bo i > 44.3). However the ASCA 
data show a strong bias at the low luminosity end of the distribution due to the absence of 
objects in the database. We believe this selection effect accounts for the apparent flattening 
of the L-T relationship for the high z clusters at low luminosities (Figure 1). Thus we 
restrict our analysis to 45.2 < log 10 L bo i < 45.7 for which there are no strong selection effects 
and for which there is maximum overlap of the low and high z data sets and a reasonable 
set of objects (39). Applying the two dimensional Kolmogorov-Smirnov test we find that at 
45.2 < log 10 L bo i < 45.7 the low and high z data show no evidence of being significantly 
different (probability of being different < 60% for all high/low redshift divisions). The 
inclusion of the lowest luminosity points, combined with a smaller data set led Tsuru et al 

(1994) to infer a change in the L-T relation which is not seen in the present data. 

The high z, large accretion rate, cooling flow clusters in Figure 1 (e.g. RXJ1347, 
A1835, EMSS1455 occupy a ridge line to the right in the L-T plot. This property was noted 
in the low z sample by Fabian et al (1994). 

We estimate the allowed change in either the normalization or slope of this relation 
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by examining the variation in the mean temperature in redshift shells. In Figure 2 the 
mean cluster temperatures determined by maximum likelihood are plotted as a function 
of redshift bin. Much of the plotted allowed range in the mean temperature is not due to 
uncertainties in the data but to real width in the L-T relation, as seen in the low z sample 
(Figure 1 and [Fabian et al 1994 ). Solid vertical error bars represent the 90% confidence 



limits on (T), dotted error bars indicate the intrinsic dispersion of the L^-T distribution 
estimated by likelihood analysis. Using go — 0.5 gives a small but significant change in the 
normalization, such that higher redshift clusters are less luminous for a given temperature. 

The limits on temperature evolution are estimated from Figure 2 as: Alog 10 (T) ~ 0.04 
for Alog 10 (l + z) ~ 0.15 at fixed luminosity. 

We also compare the relationship of velocity dispersion to temperature and luminosity, 
using the large sample of Fadda et al (1996) for the low z sample and the Carlberg et 
al (1996) and Fabbricant et al (1991) data for the high z sample. We find that the high 
z clusters in the temperature range 4-9 keV show virtually identical oy-T and oy-L 
relationships to those of the low z sample (Figure 3). AC118 and A1689 are clearly 
discrepant, indicating non-virial velocity dispersion. 



4. Discussion 

The lack of evolution in the L-T relation at z < 0.5 combined with the recent upper 
limits on the evolution of the luminosity function over the same redshift range places strong 



constraints on all models of cluster evolution (|Kaiser 1991|) . As pointed out in detail by 
Navarro, Frenk and White (1995) a model in which the intial specific entropy is large 
results in little, if any, evolution in the L-T plane and better matches the overall slope of 
the observed L-T relation. Metzler and Evrard (1994) describe a model in which entropy 



- 8- 



due to the creation of the metals is explicitly included, increasing the overall entropy of a 
kT~ 6 keV cluster by a factor of 2 and producing a value of the entropy that agrees well 
with that of several nearby clusters. This result is sensitive to the thermalization of the 
supernova shock energy. However, as David et al (1996) point out, the central entropy does 
vary from cluster to cluster, in agreement with the scatter in the L-T plot. 

Bower (1997) has parameterized the evolution of cluster central gas entropy as a power 
of the expansion factor: s min = s min (z = 0) + c v eln(l + z) (c v is specific heat capacity of 
the gas at constant volume). In this case e = corresponds to constant entropy. Data on 
the evolution of the cluster XLF and L-T relationships can then be used to constrain two 
parameters; e and n, the power spectrum index of matter density fluctuations ( [Bower 1997| ) . 
Our results imply that over the redshift range < z < 0.4, e = ± 0.9. If the evolution 
of the ICM followed that of the dark matter (e.g. [Kaiser 1986Q then our limits on e would 



imply n > 2 always, thereby firmly ruling out self-similar evolution of clusters. If there 
is little or no evolution in the XLF (e.g. |Ebeling et al 1997| ) then our result implies that 



— 1.5 & n — 0.5. Recent measurements of the low luminosity cluster population at high 



redshift QJones et al 1997| ) are also consistent with zero evolution in the XLF. 



Bryan (1996) presents a scaling relation between Lboi and T with a (l + z) 3 ^ 2 dependence 
for Q = 1. The dashed curve in Figure 2 represents this scaling, normalised to the luminosity 
band used here. Clearly the Bryan model cannot be rejected, although it does not appear 
to describe the data particularly well. The same statement can be made about the results 
of Kitayama and Suto (1996) who have improved on Press-Schechter theory by including 
the epoch of cluster formation, and predict a similar (negative) temperature evolution. 

The lack of evolution in the oy-T relation also indicates that the gas and the galaxies 
are sampling the same potential (c.f. |Carlberg et al 1996| ) and that both change in the 



same way with redshift, a result not generally expected in many models. 
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Presumably a similar lack of evolution would be seen in low density universe models 
in which cluster evolution occurs early and structure stops forming at z ~ Q^ 1 . However, 
detailed calculations for such models are not yet available. 



5. Conclusion 

We find no evolution in the cluster L-T relation as a function of redshift. These 
results are consistent with models in which the cluster core acquired a high initial entropy 
( [Navarro, Frenk and White 1995|) as required by the recent determination of the origin of 
the cluster metals flMushotzky et al 1996|) and the lack of evolution in the cluster metallicity 



( [Mushotzky fc Loewenstein 19971) . However, it is likely that this result is also consistent 



with low density cosmologies in which clusters stop evolving at relatively high redshifts, 
or other variants of the standard models. We hope that this large data set will enourage 
detailed theoretical modeling. 



We thank Y. Tanaka for early encouragement of this project. This research has made 
use of data obtained through the High Energy Astrophysics Science Archive Research 
Center Online Service, provided by the NASA/Goddard Space Flight Center. 
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Fig. 1. — X-ray luminosity vs gas temperature for low redshift (z < 0.1) (open squares) 
and high redshift (z > 0.14) (dark triangles) clusters. The errors in the gas temperature 
are symmetrized 90% confidence errors. The dashed and solid lines are unweighted linear 
regression fits to the high and low z data respectively. 

Fig. 2. — The mean cluster temperature (in the luminosity band 45.2 < log 10 L bo ; < 45.7) 
estimated by maximum likelihood is plotted versus redshift. Redshift bins contain 6 objects, 
the highest z bin at z ~ 0.45 contains 3 objects. Solid, vertical, errorbars show the 90% 
confidence limits on the mean, horizontal errorbars show size of redshift bin. Dotted, vertical, 
errorbars show the intrinsic variance in T, estimated by maximum likelihood. The dashed 
curve is the Bryan (1996) model prediction normalized to this luminosity band. 

Fig. 3. — The optical velocity dispersion vs X-ray temperature. The errors in temperature 
are as in Figure 1 while the errors in the velocity dispersion are symmetrized 68% confidence 
errors. We have labelled A1689 and AC118 as being obvious outliers of the distribution 
function. The line drawn assumed that kT=cr, 2 / and is not a fit to the data 
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Note. — Numbers in brackets in final column correspond to the 90% confidence limits on 
the temperatures. An inspection of other works (e.g. pavid ct al 1993 , Ebcling ct al 1996 ) 
indicates a scatter of £ 20% in quoted L bo i which we attribute to differences in techniques. 
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